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Abstract. While it has been shown repeatedly that ocean 
conditions exhibit an important control on the behaviour of 
grounded tidewater glaciers, modelling studies have focused 
largely on the effects of basal and surface melting. Here, a 
finite-element model of stresses near the front of a tidewater 
glacier is used to investigate the effects of frontal melting on 
calving, independently of the calving criterion used. Appli- 
cations of the stress model to idealized scenarios reveal that 
undercutting of the ice front due to frontal melting can drive 
calving at up to ten times the mean melt rate. Factors which 
cause increased frontal melt-driven calving include a strong 
thermal gradient in the ice, and a concentration of frontal 
melt at the base of the glacier. These properties are typical of 
both Arctic and Antarctic tidewater glaciers. The finding that 
frontal melt near the base is a strong driver of calving leads 
to the conclusion that water temperatures near the bed of the 
glacier are critically important to the glacier front, and thus 
the flow of the glacier. These conclusions are robust against 
changes in the basal boundary condition and the choice of 
calving criterion, as well as variations in the glacier size or 
level of crevassing. 



1 Introduction 

The calving of icebergs makes up a large component of the 
mass balance of many polar ice sheets and glaciers (Hagen 
et al., 2003; Rignot and Kanagaratnam, 2006). Recent obser- 
vations of increased calving and coincident flow acceleration 
at a number of tidewater glaciers (Joughin et al., 2004; Luck- 
man et al., 2006; Stearns and Hamilton, 2007; Howat et al., 
2007) have shown that calving and glacier flow are strongly 
linked. However, the modelling of calving processes is still 
problematic, leading to a great deal of uncertainty in predic- 
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tions of the future behaviour of tidewater glaciers, and their 
consequent contributions to sea-level rise. 

Reviewing the mechanisms of calving at tidewater 
glaciers, Benn et al. (2007b) distinguish between first-order 
calving mechanisms associated with longitudinal stretching 
of the glacier, and second-order mechanisms associated with 
buoyant forcing at the front. A calving law based on sim- 
plified physics (Benn et al., 2007a) has shown some success 
in replicating these first-order processes (Mottram and Benn, 
2009; Nick et al., 2010; Cook et al., 2012). However, no 
model exists in current practice to account for the second- 
order processes, despite the possibility that these may ac- 
count for the majority of calving in some circumstances. 

It was first noted by Weertman (1957) that there is a nec- 
essary imbalance at any ice front between the glaciostatic 
pressure outwards and the hydrostatic pressure inwards. In 
floating ice this imbalance manifests itself as a pure bending 
moment acting on the ice front, while in grounded ice it par- 
tially expresses itself as a net outward force. Reeh (1968) 
demonstrated the consequences of this effect on an analytic 
model of an ice shelf with a Newtonian rheology. He showed 
that a maximum in both tensile stress and surface elevation 
develops about one ice thickness from the front, and that this 
stress leads to calving, with the calving rate determined by 
the glacier's thickness and viscosity. He also noted the possi- 
ble effects of variations in the shape of the ice front, although 
these were not incorporated into his model. 

This same stress maximum has been described by a num- 
ber of authors. Fastook and Schmidt (1982) showed the same 
effect in a finite-element model of a Newtonian floating ice 
shelf, while Hanson and Hooke (2000) demonstrated a sim- 
ilar stress maximum in a grounded model with a non-linear 
rheology. Scambos et al. (2009) showed similar results for 
the flexure of an ice shelf, using realistic temperature pro- 
files as inputs to a Glen-type rheology. 

While the idea that frontal melting could be a driver of 
calving has been mentioned by a number of authors (Han- 
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son and Hooke, 2000; Vieli et al., 2002), there have been few 
quantitative studies of the phenomenon. A number of au- 
thors (Kirkbride and Warren, 1997; Benn et al., 2001; Rohl, 
2006) have identified the melting and erosion of waterline 
notches in lake-terminating glaciers as controls on their calv- 
ing rates. Similarly, estimates of frontal melting at LeConte 
Glacier, Alaska (Motyka et al., 2003) and several Green- 
landic glaciers (Rignot et al., 2010; Sutherland and Straneo, 
2012) show that melting may be an important term in the 
frontal mass balance of these glaciers. Jenkins (2011) has 
provided some explanation of these results, using a model of 
plume-driven melting. Given the evidence for the importance 
of frontal melt as a mass balance term, and the likelihood of a 
connection with calving processes, it has become a subject of 
interest how tidewater glaciers react to melt-driven changes 
in their front geometry. 

In this study a simple model of two-dimensional laminar 
ice flow is described, and its results are used to investigate 
the response of a tidewater glacier to submarine frontal melt- 
ing. Section (2) outlines the model equations, boundary con- 
ditions and implementation. The theoretical framework in 
which the results are presented is described in Sections (3) 
and (4), which introduce the key ideas of stress retreat and 
wet calving multiplier, and uses model results to show that 
the latter of these is well-founded. Finally, Section (5) shows 
how these can be used to predict the changes in calving be- 
haviour that should be expected under particular conditions 
of frontal melt. 



2 Ice flow model 



sors. The deviatoric stress tensor is denoted here by r^, and 
the second invariant of this tensor calculated as: 



(3) 



I ( fr 1 + 'sx 7 ) can l ^ en ^ e expressed 
as a function of the deviatoric stresses: 
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The quantities A and n are considered parameters of the 
flow law. In this work n is fixed at its conventional value of 
3, and A is allowed to vary as a function of temperature T: 



A = A cxp(-Q/RT) 



(5) 



The universal gas constant R is fixed, while the activa- 
tion energy Q and the multiplicative factor A a are chosen to 
match measurements. For the purposes of this study, stan- 
dard values are used for Q and A , as given by Paterson 
(2000). 

To translate between the stress-strain relationship of Equa- 
tion (4) and the viscosity-based Stokes formulation, Glen's 
law is inverted, giving the effective viscosity in terms of the 
strain rates and the flow parameter A: 



(6) 




2.1 Flow of ice 



2.2 Boundary conditions 



While ice exhibits visco-elastic properties on hourly to daily 
time-scales, the motion of glaciers over time-scales longer 
than about a day is best described as that of viscous fluid 
flow (Paterson, 2000). As the Reynolds numbers involved 
are extremely low, of the order of 10" 13 , and the ice density 
is constant throughout most of the glacier body, the flow can 
be approximated using the incompressible Stokes equations: 



Vp = {iV 2 u + f 
V-u = 



(1) 
(2) 



Here p is the scalar pressure, u is the flow velocity, / is 
the body force (in this case gravity), and /i is the dynamic 
viscosity of the fluid. As ice exhibits non-Newtonian flow 
properties, \x is itself a function of the velocity field, and must 
be specified through a flow law. 

The most common flow law in use for glaciological pur- 
poses is that derived by Glen (1952). This is most simply 
expressed as a relationship between the stress and strain ten- 



On the upper surface of the glacier, as well as that portion of 
the front which is above the waterline, a stress-free boundary 
condition is applied. Below the waterline, hydrostatic pres- 
sure is applied. At the rear of the domain, far from the front, 
the flow velocity is set to a constant of 1 km a" 1 . Sensitiv- 
ity tests (not shown here) indicate that the results are largely 
insensitive to this value. 

The basal boundary condition for glacier flow is an active 
topic of research, and several relationships have been pro- 
posed. The most common of these are those based on the 
work of Weertman (1957), who gives a power law relation- 
ship between basal shear stress T(,, basal velocity u\>, and the 
effective pressure at the base N = Pi—p w : 



TV 



■ku p b N- q 



(7) 



This relationship, while convenient numerically, and diffi- 
cult to disprove empirically, has been shown to have difficul- 
ties, notably the lack of an upper bound for shear stress. A 
calculation by Iken (1981) showed that for a bounded basal 
slope, such a bound must exist. This result was reproduced 
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in a more general setting by Schoof (2005), who suggested a 
relationship which was refined by Gagliardini et al. (2007): 



n 



\*AC n N n 
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(8) 



Here A* and m* are the dominant wavelength and slope 
of bed features, and C is a constant subject to the inequality 
C <m* . Following Pimentel et al. (2010) the relationship 
C = 0.84m* is used, based on the result for a sinusoidal bed. 
We assume that water pressure is hydrostatic, based on a free 
connection to the ocean. 

2.3 Numerics 

For the purposes of this study, a two-dimensional solution to 
the incompressible Stokes equations is sought, using Glen's 
flow law as a constitutive relation. The glacier is treated as 
of uniform width, and sufficiently wide that lateral boundary 
effects are unimportant. 

Such a solution is produced using the free open-source fi- 
nite element solver FreeFem++ (Hecht et al., 2005), using 
a standard triangular P2 element for the velocity field and a 
PI element for the pressure, both implemented on an unstruc- 
tured grid. A finite element implementation allows the model 
to easily handle a variety of geometries, as well as allowing 
the model's resolution to be focused on the areas of interest. 
In order to handle the implicit definition of /i through Equa- 
tion (6), the system is solved iteratively, beginning with the 
Newtonian solution. 

The basal boundary condition (Equation (8)) is also non- 
linear, which presents some difficulties in the numerical im- 
plementation. A Robin-type boundary condition is used, ex- 
pressing Tj, as a linear multiple of u at the base, and recal- 
culating the constant of proportionality with each step of the 
non-linear iteration procedure. This iterative process is com- 
bined with that for the effective viscosity, in order to reduce 
the total number of iterations required. 

This procedure is found to be much more numerically 
stable than the alternatives, such as Dirichlet or Neumann 
boundary conditions, and it does not substantially increase 
the number of iterations required for convergence over the 
case of a fixed basal velocity. 

3 Stress retreat 
3.1 Assumptions 

The aim of this study is to quantify the effects of submarine 
frontal melting on calving rate. To this end, some assump- 
tions must be made about the calving criterion. It is assumed 
that calving behaviour is determined solely by the viscous 
properties of ice, neglecting any effects due to elastic defor- 
mation. This greatly simplifies the calculation of stresses, 
and is in line with most glaciological practice. 



It is also assumed that the variable of greatest interest is the 
(Cauchy) stress, rather than the strain rate (or equivalently, 
deviatoric stress). This is supported by Vaughan (1993), who 
notes that strain rates at crevasse sites vary by almost three 
orders of magnitude while stresses are almost constant. For 
simplicity, only the first (most tensile) principal stress is con- 
sidered, as it is assumed to be the controlling factor on frac- 
ture. External sources of stress are ignored, based on the 
analysis of Bassis et al. (2008). 

It is known (Weertman, 1973) that water-filled crevasses 
are likely to penetrate the full thickness of the glacier, un- 
der most realistic stress conditions. Therefore, if we assume 
that crevasses which reach below the waterline are likely to 
contain water, we can disregard the stresses which occur at 
depths much below the waterline, for the purposes of calcu- 
lating crevasse depths. 

Similarly, if we assume that the dominant control on the 
crevasse's growth is the stress field around the crack tip, the 
stress field near to the surface of the glacier is irrelevant, once 
the crack has grown beyond a certain size. Thus, assuming 
that compressive stress increases with depth, the most impor- 
tant factor in determining whether the crevasse penetrates the 
whole glacier is the stress field around the waterline. 

It should be noted that the Benn-Nye calving criterion 
(Benn et al., 2007a), as widely used in models of tidewater 
glaciers, fits all of the above assumptions, as do any obvious 
modifications of it. 

Finally, it is assumed that internal deformation of the 
glacier takes place on timescales much longer than those as- 
sociated with calving or frontal melting. This follows from 
the calculation of van der Veen (2002), who showed that a 
typical tidewater glacier front deforms at a rate which is be- 
tween one and two orders of magnitude too slow to be re- 
sponsible for calving. 

It therefore suffices to look only at the cumulative amount 
of frontal melt, rather than the rate at which it occurs. An- 
other advantage of this assumption is that it can be internally 
verified by using the stress field to calculate instantaneous 
velocities. Note that no assumption is made about the abso- 
lute velocity of the glacier, merely about the rate at which it 
is deforming internally. In fact, for a stable calving front it is 
necessary that the glacier sliding velocity is of a similar or- 
der to the calving rate, and thus must take place on the same 
timescale. 

3.2 Definitions 

In this study, the dimensionless dry calving length e is de- 
fined as the aspect ratio of the incipient iceberg, when no 
frontal melt is occurring. Thus, for a glacier of thickness H, 
the distance from the ice front to the point of calving is eH. 
While this measure is undoubtedly variable, and may even 
include a stochastic component, it is fixed for a given calving 
event. It is also useful to define the undercut length d as the 
depth-averaged cumulative frontal melt, and the dimension- 
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Fig. 1. Above: Schematic of the geometry used in this study. Be- 
low: Schematic of change in stress due to an undercut. The stress 
retreat r is the distance between a stress contour in the baseline 
state, and the equivalent contour in the perturbed state, measured 
along the waterline. 



those same conditions exist at the point x + r = x+ujd. Thus, 
if calving is occurring at an interval At, leading to a dry calv- 
ing rate of a;/ At, the calving rate once frontal melting is in- 
corporated will be equal to xj At + cod/ At, an increase of w 
times the mean melt rate. 

It may therefore be useful to treat lj as a measure of the 
sensitivity of the calving rate to variations in frontal melt rate, 
with a value of one corresponding to a simple additive model. 
Such a way of thinking about the effects of submarine melt 
is the dominant one in those works which include melting 
in the frontal mass balance calculation (Motyka et al., 2003; 
Amundson and Truffer, 2010). While such an approach is 
observationally correct — there is no easy way to distinguish 
frontal melt-driven calving from any other kind — it will 
be found more productive here to make this distinction, and 
to separate 'dry' melt-free calving from 'wet' frontal melt- 
driven calving. 

This interpretation must be used with care, however. An 
increase in the mean size of a calving event is likely to have 
repercussions on the dynamics of the glacier. As such, the 
effects described here should be considered merely one com- 
ponent of a system of interacting processes and feedbacks 
which ultimately determine the behaviour of the glacier. 



4 Wet calving multiplier 



less waterline height h as the ratio of the glacier depth below 
sea level to the full glacier thickness. 

The central concept of this analysis is the 'stress retreat', 
which is a measure of the spatial effect of a stress perturba- 
tion, defined as follows. A reference frame is used whereby x 
is a horizontal variable, increasing inland from zero at the ice 
front. Given a reference stress state a re f(x) and a perturbed 
state a pert {x), the stress retreat r is the minimum distance 
inland such that a re f(x) = a pert (x + r) . In other words, the 
effect of the perturbation is to move the stress field inland by 
a distance r. This distance is, of course, a function of posi- 
tion x, and may vary considerably. In all cases considered 
here, r is positive and finite. 

Intuitively, the stress retreat can be thought of as the dis- 
tance that the stress field is 'pushed back' by a perturbation, 
such as undercutting or a change in the frontal boundary con- 
dition. Ice in this situation will behave as though it were this 
distance further forward in an unperturbed glacier. 

Finally, the wet calving multiplier to is defined as r/d, 
where r is the stress retreat due to an undercut length d. From 
first principles, there is no reason to suppose that to is inde- 
pendent of d, but this shall be shown empirically to be the 
case in Section (4). 

The ratio u> can be interpreted as follows: assuming, in 
the absence of frontal melting, that the conditions for calving 
exist at a point x, then after a quantity of frontal melting d, 



In order to calculate the effects of undercutting on near- 
frontal stress, the flow model is used to compute stress fields 
in a variety of configurations. As a baseline unperturbed 
model run, meant to simulate a typical medium-sized tidewa- 
ter glacier, the model is run in a flat rectangular slab configu- 
ration. The ice thickness is 300 m, and the basal parameters 
are A* = 20 m and m* = 0.13 — see below for the sensi- 
tivities to these parameters. The water level varies between 
model runs, as the results are quite sensitive to this variable. 

By altering the shape of the domain, the effect of undercut- 
ting by frontal melt can be simulated. Figure (2) shows the 
deviatoric stress fields generated by a selection of geome- 
tries and water levels, assuming uniform frontal melt below 
sea level, and an isothermal glacier. Qualitatively, it seems 
clear that undercutting results in an increase in tension due 
to the bending moment exerted by the overhang, as well as 
the reduction in basal traction near the glacier foot. These ef- 
fects increase with the undercut length d and show qualitative 
variation with the water level h. 

Figure (3) shows the stress retreat as a function of under- 
cut length for a variety of scenarios. In each case, a uniform 
undercut is introduced and the stress retreat measured, rela- 
tive to a particular dry calving length e. For each scenario, 
a linear fit through the origin is possible, with R 2 > 0.99. 
The slope of this fit is equal to the wet calving multiplier to, 
which is henceforth assumed to be independent of the under- 
cut length. Although this assumption must certainly break 
down at large undercut lengths, it appears to hold for under- 
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Fig. 2. Comparison of first principal Cauchy stresses, relative to 
hydrostatic pressure, for a variety of water depths and undercut 
lengths, assuming uniform undercutting. Grey dashed line indicates 
water level. 
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Fig. 3. Stress retreat as a function of undercut length for a variety of 
scenarios involving different water depths h and dry calving lengths 
e. In all cases the relationship is very close to linear. The slope of 
the linear fit is lj, the wet calving multiplier. 



cut lengths which are up to 20% of an ice thickness, far larger 
than the expected depth of real- world undercuts. 



5 Sensitivities 

Given that lj is well-defined for a given scenario, the ques- 
tion remains as to what factors influence its value. The most 
obvious of these is the dry calving length e. As this is used 
as the initial point from which the stress retreat is measured, 
it should come as no surprise that the magnitude of the stress 
retreat (and hence uS) is dependent on its value. 

For a grounded or partially grounded glacier, values of e 
greater than one are usually held to be unlikely, given that the 
resulting berg would be unable to capsize, and would thus 
have no obvious route of escape from the glacier. Here, the 
upper limit is drawn at a value of e = 1.5, to allow for some 
leeway in the system. Similarly, values of e < 0.25 are ne- 
glected, as such narrow calving events are likely to be much 
more affected by the detailed geometry of the front than by 
viscous stresses. However, it should be noted that they are 
not ruled out by this model, merely likely to be modelled 
incorrectly. 

Another variable of interest is the water depth, or more 
loosely, the 'degree of grounding'. As there is known to be a 
significant difference in calving behaviour between grounded 
and floating glaciers (Walter et al., 2010), it seems reasonable 
to suggest that the water level may have a significant qualita- 
tive effect on calving, even if the transition is not as abrupt as 
that between grounded and floating ice. As such, a selection 
of water depths are investigated, ranging from h — 0.5, for 
a well-grounded glacier, to h = 0.85, a glacier almost at the 
point of flotation. Flotation occurs at h = Pi/p w — 0.89, at 
which point the calculations of the model diverge in a non- 
useful manner. 

As can be seen in Figure (4), typical values of lj for the 
simplest scenarios are in the range 1—4, indicating that in this 
idealised situation, frontal melt drives calving at a rate up to 
four times the mean melt rate. Higher water levels generally 
lead to larger values of the multiplier lj, as do shorter calving 
lengths. In the h = 0.5 case, where the glacier is immersed in 
water to its midpoint, the multiplier effect is relatively weak, 
with frontal melt-driven calving occurring at around one and 
a half times the melt rate. The effect is much stronger in the 
more typical h = 0.7 case, with a multiplier in the range of 
two to three. 

The highest values of lj are usually found at e ~ 1 — h, 
indicating that the geometry of the above-water portion of 
the nascent iceberg is important. At calving lengths shorter 
than this, there is some drop-off in the value of lj, although 
as previously stated, the model may not be fully capturing 
the complexities of the stress field so close to the front. 

Glaciers close to flotation also show an increase in uu for 
calving lengths of around an ice thickness, foreshadowing 
the divergence in lj when the glacier comes afloat. However 
care should be taken in interpreting this result, as the stress 
distribution becomes very uniform here, and in reality stress 
variations from other sources will probably play more of a 
role. 
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Dry calving length (e) 

Fig. 4. Relationship of wet calving multiplier oj to calving length e 
for varying water depths. Also shown is the floating case, where oj 
rapidly diverges. 



Significantly, in all but a few cases, cu is greater than one. 
This means that undercutting could be driving calving at rates 
greater than the frontal melt rate itself, providing an ampli- 
fication of the oceanographic forcing on the glacier. This 
effect will be most strongly felt on those glaciers which are 
themselves well immersed in water. This may go some way 
towards explaining the empirical relationships which have 
been identified between calving and water depth (Sikonia, 
1982; van der Veen, 1996). 

5.1 Thermal regime and viscosity 

Given that many ice sheets are far from isothermal (Paterson, 
2000), it might be assumed that the thermal regime of their 
outlet glaciers is similarly heterogeneous. Although mech- 
anisms have been suggested by which the hydrological sys- 
tem of a glacier could result in near-isothermal conditions 
(Phillips et al., 2010), this has not been widely observed in 
the field. As such, it is useful to consider the effects of a 
thermal gradient on the frontal melt-calving relationship. 

For simplicity, a linear temperature trend within the ice 
is assumed, ranging from -10 °C at the surface to 0°C at 
the bed. The effect of pressure on the melting point is ne- 
glected for simplicity. The effects are shown in Figure (5). 
In general, for glaciers close to flotation, the effect of the 
temperature gradient is to boost oj, in some cases by around 
a factor of two. For more heavily grounded glaciers, such as 
the h — 0.5 case, the effect is very slightly negative. 



8.25 0.50 0.75 1.00 1.25 1.50 
Dry calving length (e) 

Fig. 5. The effect of a temperature gradient on the ui-e relationship. 
Solid lines show the baseline cases, while dashed lines show the re- 
lationship when a 10 °C temperature gradient is applied. For lightly 
grounded glaciers, there is a positive effect, while heavily grounded 
glaciers see almost no difference. 



The implication here is that glaciers in colder regions, with 
less developed hydrological systems and thus less homoge- 
neous temperature profiles, are likely to be more severely im- 
pacted by a constant degree of undercutting, although these 
glaciers are also likely to be in colder marine environments, 
which will tend to decrease the quantity of submarine melt. 
In fact, higher atmospheric temperatures and the consequent 
surface melt may act to stabilize glacier fronts, by reducing 
the effects of undercutting on calving, if surface meltwater is 
sufficient to render the glacier temperate near the front. 

It is likely that a two-part classification scheme is neces- 
sary here, distinguishing both between polar and temperate 
glaciers as is usual, but also between those terminating in 
warm and cold water. While there is certainly a strong cor- 
relation between the two groups, the relationship is imper- 
fect, and there are certainly examples in Greenland of 'cold' 
glaciers coming in contact with relatively warm water (Rig- 
not et al., 2010; Sutherland and Straneo, 2012; Christoffersen 
et al., accepted). It has also been shown (Seale et al., 2011) 
that the behaviour of tidewater glaciers in Greenland respects 
strong geographical boundaries, which are more associated 
with oceanic than atmospheric or glacial temperatures. 

Up until now, the glacier has been treated as a completely 
solid block of ice. In reality, the upper surface of a tide- 
water glacier is often heavily crevassed, resulting in reduced 
strength in the upper layers. Crevassing such as this is often 
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Fig. 6. The effect of surface damage on the ui-e relationship. 
Dashed lines show the relationship when an enhancement factor is 
applied to the upper layers. No discernible pattern is visible, and 
the difference is likely to be a numerical artifact. 



Fig. 7. The effect of a change in basal boundary condition on the 
uj-e relationship. Dashed lines show the relationship with an ex- 
tremely simple linear basal boundary condition. The most signifi- 
cant change is the reduction in oj for the h = 0.85 case. 



represented in flow models through the use of an enhance- 
ment factor, which reduces the effective viscosity of the ice 
in such areas. Now, an enhancement factor is applied, reduc- 
ing the effective viscosity above the waterline. At the surface 
the viscosity is reduced by a factor of ten, and this reduction 
is scaled linearly with depth until the viscosity reaches its 
normal value, 5m above the waterline. This height is cho- 
sen so as to minimise the effects of a 'kink' in the effective 
viscosity at the waterline, which interferes with the accurate 
calculation of stresses. 

The effects of this change can be seen in Figure (6). While 
there are some small differences, it is quite likely that these 
are due to numerical inaccuracies, due in large part to the 
'kink' in the effective viscosity near the waterline. In general, 
it seems that the rheology of the above-water portion of the 
glacier has very little effect on the value of uj. This is a very 
helpful result, as it provides some validation for the approach 
of using an idealized slab glacier model. 

5.2 Basal boundary condition 

The basal boundary condition is one of the most uncertain 
elements in any glacier model. Therefore, it seems prudent 
to check that the results given here are robust against changes 
in this area. As a simple test, the boundary condition given 
in Equation (8) is replaced with a simple linear relationship 
between sliding velocity and basal traction, with the constant 



of proportionality (50 Pa a m" 1 ) chosen to best match the ve- 
locities of the more complex model. 

The results are shown in Figure (7). The largest effect is 
the reduction in the value of to for the nearly floating case, 
as the changes in effective pressure at the bed no longer can 
have any effect. In the other cases, there is a small decrease 
in the value of w, but the general pattern remains very similar. 
Therefore it can be concluded that the basal boundary condi- 
tion has little effect on this particular calving mechanism. 

Another possible source of variability is the glacier thick- 
ness itself. Although it can be argued through dimensional 
analysis that most of the physics should be unchanged by 
an increase or decrease in physical scale, there are some ef- 
fects surrounding the basal boundary condition which do not 
scale in a simple way. However, sensitivity tests (not shown) 
indicate that the glacier thickness is unimportant, with the re- 
sulting graphs being visually indistinguishable. As such, the 
results given here can be thought of as independent of glacier 
thickness. 

5.3 Undercut shape 

Another consideration which has to this point been ignored is 
the vertical distribution of melt over the glacier face. As es- 
timates of frontal melt rates have generally arisen from heat 
balance calculations, there is no empirical evidence to sug- 
gest a particular form for the frontal melt profile. However, 
given the modelling results of Jenkins (2011), it should be 
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Fig. 8. Above: The u-e relationship for a variety of frontal melt 
profiles. Below: Cauchy stress fields (relative to hydrostatic pres- 
sure) associated with frontal melt profiles. 



expected that the frontal melt rate would be vertically inho- 
mogeneous. 

An attempt can be made to quantify the effect of this in- 
homogeneity by using a variety of different idealized frontal 
melt profiles. As well as the uniform profile used thus far, the 
model can be run using undercuts in wedge shapes, as well 
as a parabolic curve which peaks in the center of the subma- 
rine ice face. Note that in all cases, the undercut length used 
in calculations is the mean melt rate on the ice face. This 
is equivalent to half the maximum frontal melt rate for the 
wedge shapes, or two thirds that in the parabolic case. 

Figure (8) shows the results of these tests. The differ- 
ences are large. The basal wedge results in values of cj which 
are consistently about 50% greater than those in the uniform 
case. The surface wedge, by contrast, results in values which 
barely get above one, meaning that in this case, the simple 




0.50 0.75 1.00 1.25 
Dry calving length (e) 



50 



Fig. 9. The effect of a change in surface slope on the u-e relation- 
ship. Dashed lines show the relationship with a surface slope of 5 °. 
Note that in this case, h refers to the waterline height at the ice front 
only, with effectively smaller values of h further from the front 



additive approach is sufficient. Interestingly, the parabolic 
profile has a similar value of w to the uniform case for small 
values of e, but it drops off much more quickly for larger 
values. 

Looking at the stress distributions, it can be suggested that 
the reason for the large value of uj in the basal wedge case 
is the movement of the 'fulcrum' about which the glacier is 
bending to the base of the glacier. This area is visible in 
Figure (8) as a region of high deviatoric stress. While the 
surface wedge case includes much higher stresses, they occur 
at the waterline and thus have only local effects. 

5.4 Surface slope 

A related issue to that of the undercut geometry is the ge- 
ometry of the above-water portion of the ice, or equivalently 
the surface slope of the glacier. For simplicity of modelling, 
and due to the extremely complex and variable nature of real 
glacier geometries, we restrict our attention here to constant 
surface slopes. 

Figure (9) shows the effect of a 5 ° surface slope on the 
calving multiplier. As might be expected, this has the effect 
of reducing the multiplier, relative to a glacier with a flat sur- 
face, as the sloped glacier is effectively 'more grounded' in- 
ward of the front, a situation which we have shown results in 
a reduced multiplier. However, the effect is still significant, 
particularly at short calving lengths. 
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In the event of a reverse surface slope, as observed at 
Helheim and Kangerdlugssuaq by Joughin et al. (2008), we 
should expect a similarly reversed effect on the multiplier. 
This is likely to lead to enhanced calving at longer lengths, 
as was observed by that study. 

Owing to numerical simplifications in our model, we are 
unable to investigate the effects of a bed slope. While we be- 
lieve these effects are likely to be quantitatively significant, 
we expect that they will not change the main qualitative con- 
clusions of this study. Future work, applying this approach to 
real-world glacier conditions, will likely need to incorporate 
such effects. 



6 Conclusions 

The model described here is a novel ultra-high-resolution 
finite-element two-dimensional flowline model of the 
stresses near an ice front, incorporating state-of-the-art 
boundary conditions and full-Stokes ice physics. It provides 
a rapid and flexible means to investigate the stress fields near 
the calving fronts of tidewater glaciers. Using a simple set 
of assumptions about the form of a calving law, a theoretical 
framework has been developed which allows for conclusions 
to be drawn about the effect of a frontal perturbation on calv- 
ing. It is shown that the effects of submarine frontal melting 
can be described within this framework, and that the impacts 
on calving may be much greater than previously assumed. 

The model demonstrates that with uniform submarine 
melt, extra calving is generated at a rate between one and 
four times the melt rate. This effect is greater the closer a 
glacier is to flotation, and can be exacerbated by the existence 
of a thermal gradient within the glacier, or melting which is 
focused near the base of the ice front. Given the results of 
Jenkins (2011), which show that this form of frontal melting 
is quite likely on many glaciers, this result is of particular 
importance. 

These results provide more evidence that the effects of wa- 
ter temperature near the base of tidewater glaciers are critical 
to the behaviour of the ice fronts. These findings are robust 
against both the basal boundary condition and the choice of 
calving law, and appear to be independent of the glacier size, 
as well as the level of damage present at the glacier surface. 
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